Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 February 2008 (MN M£X style file v2.2) 



Radiative Transfer Along Rays in Curved Spacetimes 



Avery E. Broderick* 

Institute for Theory and Computation, Harvard-Smithsonian Center for Astrophysics, 60 Garden St., MS 51, Cambridge, MA 02138, USA 



4 February 2008 



in 
o 
o 

(N 
> 

o 



> 
in 

in 



in 
o 

Oh 

6 



ABSTRACT 

Radiative transfer in curved spacetimes has become increasingly important to under- 
standing high-energy astrophysical phenomena and testing general relativity in the 
strong field limit. The equations of radiative transfer are physically equivalent to the 
Boltzmann equation, where the latter has the virtue of being covariant. We show that 
by a judicious choice of the basis of the phase space, it is generally possible to make 
the momentum derivatives in the Boltzmann equation vanish along an arbitrary (in- 
cluding nongeodesic) path, thus reducing the problem of radiative transfer along a ray 
to a path integral in coordinate space. 

Key words: radiative transfer, relativity, black hole physics 



1 INTRODUCTION 



Radiative transfer in general relativistic environments is 
required to understand many high-energy astrophysical 
phenomena. Despite this, the approaches utilised are ei- 
ther ad hoc, use c omplicated transfer functions (s ee, e.g., 
Curmineha ml ll975lh moment formalisms (see, e.g., iThornd 



1981), Monte-Carlo methods (see, e.g.. lJanka fc Hillebrandt 
1989), or the Boltzmann equa t ion is solved in its en- 
tirety (see, e.g., iLindauistl Il966t iLiebendorfer et all l2005t 
and references therein). Many researchers who are inter- 
ested in producing images or light curves of compact ob- 
jects have tacitly assumed that the problem of radia- 
tive transfer is reduced t o performing an integral along 
the line-of-sight (see, e . g.. [ Jaxoszynski fc Kurpiewskilll997t 



iBroderick fc Blandforcfl 120041) . While being intuitively ob- 
vious (radiative transfer is a local process), we have been 
unable to find an explicit justification for this in the liter- 
ature. Here we show that by Fermi-Walker propagating a 
basis along the ray the Boltzmann equation may indeed be 
reduced to a path integral along the ray (as long as scatter- 
ing is ignored). 



2 BOLTZMANN EQUATION 

The Boltzmann equation is given by 



dx" dN 
dr dx^ 



+ 



dp M dN 
dr dp^ 



= S(x»,p»,N), (1) 



where r is an affine parameterisation of the ray, p M is the 
tangent to the ray, N is the distribution function of the 
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particl es under consid eration, and <S contains the source 
terms jLindauisdf igGe) 1 . The particular f orm of S a n d that 
it i s a Lorentz scalar are discus sed in iLindaui ^ *|l9^)^ 
and | Broderick fc Blandfor As shown in lLhidcmrsg 

( 1966), equation Q can be written in terms of an arbitrary 
tangent-space basis ({e^}) as 



dx" ON 
dr dxv- 



+ 



dpadN 
dr dp a 



= S(x»,p a ,N) : 



(2) 



where p a = p M e^. The momentum terms will generally 
vanish if dp a /dr = 0, which is true by definition if {e^} 
are Fermi-Walker transported 2 along the ray (see, e.g., 
iMisner et aflll973t §6.5), resulting in 



dx" dN 
dr dx^ 



= S(x»,p a ,N). 



(3) 



If the path is a geodesic, Fermi- Walker transport simply 
reduces to parallel transport. It is worth noting that the 
condition dp a /dr = will also be satisfied if is a Killing 
vector. 



3 EXAMPLE APPLICATIONS 

The formulation of radiative transfer as a path integral 
along the ray in coordinate space (as opposed to in the en- 

1 In general, it is also possible to write this in terms of the conju- 
gate coordinate-momentum pair and p M , with identical results. 

2 Fermi- Walker transport is defined by 

u'Vvf = (uV-»V)/„, 



where is the unit tangent vector [u v u v 
a M = m"V„u'\ If the ray is a geodesic a 11 ■■ 
reduces to parallel transport. 



: — 1) of the ray and 
and this trivially 
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tire phase space) is especially convenient, e.g., for modeling 
high-resol ution images and light curve s of compact objects 
(see, e.g.. iBroderick fc Loebll2005allbl(J) and ray-casting ra- 
diative transfer cod es (see, e.g.. iRazoumov fc Scottl fl999t 
ISokasian et ahlfeoOlT) . In the former, typically a congruence 
of rays are traced from a distant observer backwards in time 
towards an emitting region. In the latter an isotropic dis- 
tribution of rays are traced outward from each grid point. 
Using the radiative transfer equation in the form of equation 
generally requires that the complete distribution func- 
tion be constructed at each point along each ray. In contrast, 
equation is trivially integrated along the ray at a single 
observing frequency (which for imaging purposes is usually 
the situation of interest). This is true even when the rays 
are not geodesies, e.g., in the presenc e of strong refraction 
(see, e.g. IBroderick fc Blandfordl l20031, 

For vacuum photon propagation in the Kerr spacetime 
(that which is most likely to be astrophysically relevant), 
and in Petrov type-D spacetimes more generally, the trans- 
ported basis {e^} may be constructed algebraically. This is 
done by first noting that the tangent vector to the ray is al- 
ready parallely propagated by definition, and thus choose 
6q tx da; M /dr. Two fiducial vectors, and e%, orthogo- 
nal to each other and eg, are then chosen at some position 
along the ray. These may be transported along the ray alge- 
braically by using the orthogonality relations, a normaliza- 
tion condition, and the c omplex const ant first described by 
IWalker fc Penrose! jl97fJ) (though see IChandrasekharll992l 
as well) 3 . Finally, the fourth basis- vector is then obtained 
from the conditions 



= 1 . 



= 0. 



(•») 



If desired, an orthonormal basis can be created from {e^} 
by making a linear combinations of eff and e%. 

However, for many applications, the explicit construc- 
tion of {e^} is unnecessary. For example, consider the case 
when the source terms S depend only upon JV, position 
(x M (r)) and the momentum of photons moving along the 
ray (p m (t) oc dx^/dr), i.e., anisotropic scattering into and 
out of the ray may be neglected. Thus, only vector e(J , as de- 
fined above, is required. However, this is presumably already 
available from the ray construction. 

An astrophysically example in which this is the case is 
self-absorbed synchrotron emission. In this case S = j — ajV 
where for a power-law electron distribution, 



. n e / u>b sin i 
3 = Jo— 

LU \ LO 



a = aon e 



lu b sin ■& 



a+3/2 



(6) 



where jo and ao are constants, a is the spectral index, n e 
is the proper electron number density, lo and lob are the 



Figure 1. The intensity map of a stationary, non-thermal, 
self-absorbed synchrotron emitting plasma surrounding a non- 
rotating black hole as seen at infinity for a single frequency. The 
contrast is linear ranging from maximum intensity to vanish- 
ing intensity (black). The density (oc r -3 / 2 ) and magnetic field 
strength (<x r -5 / 4 ) are spherically symmetric and the latter is as- 
sumed to be randomly oriented (i.e., the emission and absorption 
have been averaged over field orientations). The spectral index of 
the emitting electrons is taken to be 2. For reference, the white 
ring shows the size of the horizon. 



photon and cyclotron frequency in the plasma rest frame, 
respectively, and $ is the angle between the magnetic field 
and the ray in the plasma rest frame. Clearly, along the ray 
n e and ujb are functions of x m (t) only. The photon frequency 
as measured in the rest frame is given by (where u M is 

the plasma four- velocity) and thus is a function of ^(t) 
and p M (r) only. Finally, the angle 1? may be defined by 
its cosine, which is in turn given by p ,J, u v *F l _ lv /luljb where 
*F^ V is the dual of the electromagnetic field tensor. Thus, 
as claimed S is a function soley of jV(t), x^(t) and p M (r). 
This may then be inserted into equation which can now 
be trivially integrated to yield JV(pq), where p£ enters via 
the initial conditions of the ray (presumably at t = 0). 
An image produced by this procedure is shown in Figure 
and additional more complex e xamples can be found 
elsewhere in the literature (see, e.g.. IBrode rick fc Blandfordl 
I 2003L l2004l: IBroderick fc Loedl2005 JbldlFalcke et alJl200fI 
Jaroszvnski & Kurpicwski 1997). 



3 For the Kerr spacetime, in Boyer-Lindquist coordinates the 
Penrose- Walker constant is given by (e.g., for ef) 

Kpw = 

(r — ia cos 9) j (e^e^ — egej) + a (^egef — e\\ sin 2 9 
-i [(r 2 + a 2 ) (e M - ege?) - a (e^e? - egej)] sin^} . (4) 



4 DISCUSSION 

The simplification of the Boltzmann equation has a straight- 
forward physical interpretation. Since radiative transfer is a 
local process, curvature enters only in relating the tangent 
space at different points of the spacetime. For the observer 
moving along a ray, by definition, this is naturally accounted 
for by Fermi- Walker transporting a tangent space basis, thus 
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reducing the problem to its flat space analogue (for which 
the momentum derivative terms can always be made to van- 
ish). 

While it is possible to define the basis {e^ } on any three 
dimensional hypersurface of the spacetime, unfortunately, it 
is not generally possible to produce a well defined basis in 
the full spacetime as a result of ray intersections. This does 
not mean, of course, that it is not possible to integrate the 
Boltzmann equation along any single ray. Rather, the phase 
space bases for different rays which intersect will in general 
produce a different basis at the point of crossing, and thus 
cannot be compared directly (though they may be rotated). 
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